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T he c ondensate fraction of a homogeneous and dilute Bose gas is expanded as a power series 
of V na 3 as No/N = 1 — ci(na 3 ) 1/ ^ 2 — C2(na 3 ) — c^(no?) 3 ^ 2 .... The coefficient c\ is well-known as 
8/3-i/tt, but the others are unknown yet. Considering two-body contact interactions and applying 
a canonical transformation method twice we developed the method to obtain the higher order 
coefficients analytically. An iteration method is applied to make up a cutoff in a fluctuation term. 
The coefficients ares c 2 = 2(ty - 8/3) and c 3 = (4/V?) (tt - 8/3) (10/3 - tt). 
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I. INTRODUCTION 

Homogenous and dilute (HD) Bose gas has been studied for more than a half century since the classical paper by 
Bogoliubov flj. Despite its simplicity and significant published research, there still remains some unsolved fundamental 
problems. The condensate fraction, or the ratio of particles of zero momentum in the perturbed ground state, is a very 
fundamental concept of a many-body Bose system. At zero temperature, the properties of HD gases are expressed in 
terms of a single parameter, the s-wavc scattering length a, and a quantum loop expansion is expressed by dilute gas 
parameter na 3 <C 1 where n — N/V is the particle number density. 

In a HD system the condensate fraction is expressed as @, Q 

^ = 1 - cxtna 3 ) 1 / 2 - c 2 (na 3 ) - c 3 (na 3 ) 3 / 2 - . . . . (1) 

The coefficient c\ is well-known as ci = but the others are unknown yet. In this paper, we will calculate the 

unknown coefficients Ci and C3 of the next higher order terms analytically using a double canonical transformation 
and an iteration method. 

II. BOGOLIUBOV HAMILTONIAN 

If the interparticle interaction is a contact potential, the effective Hamiltonian for the identical particles of mass m 
is written as the Bogoliubov form 

H = ^e k a\a k + — ^ al^a^akj^+k^ks+ki, (2) 

k k i k 2 

k 3 k 4. 

where e k = h 2 k 2 /2m and g is the coupling constant given by g = 4:Trh 2 a/m. a ki and a) k satisfy the bosonic commutator 

relation: [ak, atr] — o"k,k' ■ Note that this Hamiltonian is valid only for the purpose of applying first-order perturbation 
theory. 

Since No is a very large number, the operators ao and a can be regarded as \/Nq. In a two-body interaction, the 
interacting part of the Hamiltonian in Eq. © is expanded to keep the terms up to the order of N% , No, and y/No as 

aa@ 

£r 9 f t 

lint — 7^7 a a 0°0 a 



+ 



k^a 

V ^ [ a k+ q a o a ka q + a{ +q al q a k a . (3) 



k,q^0 
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The three-body interaction is not considered. The last term, which is generally neglected in textbooks, originates 
from the interactions of particles out of and into the condensate and the key of the calculation. 
Let's introduce a new variable as 

lk = ^2a\ +q a q . (4) 

Since the Nq is expressed by the number relation 

N ° = & - \Y.( a i ak + a -k a -k), (5) 

we can rewrite the model Hamiltonian as a function of the density n as 

H = \n 2 Vg 
2 J 



2 [( £fc + n 9)( a l a k + a- k a -k) + n 9( a W-k + a ka-k) 



MO 

+ 7=E^(^+ a -fc)- (6) 
ViV mo 

III. DOUBLE CANONICAL TRANSFORMATION 

Next, we apply the conventional Bogoliubov transformation of the new the operators a k and a\ to diagonalize the 
Hamiltonian p|. 

Gk = /T^2 ^ bk + Akb -k)' 

a k = it A A b i+ A kb-k). (7) 

It is clear that b k and b\ satisfy the same commutator relations as a k and at- Then, the Hamiltonian in Eq. © 
becomes 

H = ln 2 Vg + J2^T2 [(e k + ng)A 2 k +ngA k ] 
1 mo k 

+ o J2 T^J2 [i £ k + ng){l + A 2 ) + 2ngA k ] (b\b k + b f _ k b^ k ) 
1 MO k 

+ \ H 7^72 Wk + ng)A k + ng(l + A 2 k )] {b\bl k + b k b_ k ) 

MO k 

A k is chosen to make the off-diagonal term vanish. 



2(e k + ng)A k +ng{l + A 2 ) = Q, (9) 

and we obtain A k (= A- k ) as 

A k = ^-fo+m) (10) 
ng 

where 

E k = ^/{e k +ngf-{ngf. (11) 
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The E k is the Bogoliubov form of the dispersion relation. 

Substituting A k into the Hamiltonian in Eq. l[8]). we obtain the following compact form of the Hamiltonian 



fc^O 



where 



G k = 



ng (1 + A k )"f k 



N v^Af 



(12) 



(13) 



The first three terms in the right hand side of Eq. (|12p are already known and the last one is the new contribution 
beyond the textbook results 0, H, 0]- We can regard the system as a collection of quantum mechanical harmonic 
oscillators exposed to additional forces given by linear terms. The linear terms b k and b\ can be made to eliminate 
by another linear transformation of the form 

b k = c k + a k , 

b\ = 4 + a_*, (14) 

where the new operators c k and ct satisfy the bosonic commutator relations. Note that a k commutes with c k and 

a k = ct- k since b k and b\ are bosonic operators. 

Substituting the new operators in Eq. fTJJ into Eq. (fT2"l) . we can rewrite the Hamiltonian in Eq. (fT2|) as 



H = \n 2 Vg + + n 9) A l + n dM 



( E kc f k c k + E k a- k a k + 2G k 



a k 



fc#o 

We choose a k to make the linear terms vanish as 



\{E k a- k + G k )c k + (E k a k + G- k )c\ 



ak 



G- k 
E k 



(15) 



(16) 



Then, the Hamiltonian in Eq. (|15[) is now diagonalized in terms of the new operators as follows 

G k G- k 



H = -n 2 Vg + Y YZ ~£ ^ £k + n ^)4 + ngM] 

k^Q k 

= H g + Y E knk- 

k^O 



fc#0 



E k 



Hg creates the well-known ground state energy density [1, H, H[ but we will not discuss it here. 



Eg 

V 



128 ( na 
15" 



(17) 



(18) 



IV. CONDENSATE TRACTION 



The condensate fraction is defined as (N — J2k^o n k)/N . We obtain the particle distribution n k of the dilute system 
from Eqs. and (JT4J) as 

n k = {0\a\a k \Q) 

- * ■ < 1 + ^ 2 <>_ m |o) 



i-4 i-4 



1 ( e k + nq \ 1 n 2 q 2 ei , , . , , 

,f^-l) + ^^(0|7 fe 7- fc |0). (19) 
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Note that (1 + A/-)/(l — A k ) = e k /E k . The magnitude of nge k /E^, is order of 1. The first term is shown in the 
textbooks and produces the well-known c\ in Eq. |T]). On the other hand, the relation between the second term and 
n k will create c%. 

The dominant part of the unknown quantity (0|7fc7_fc|0) is obtained in the following way. Applying the definition 
of 7/j in Eq. ([4]), the ground state average is written as 



(0|7fc7-fc|0) 



i E 

q,q' #o 



4+ f A a -fc+ g ' a 9'l >- 



(20) 



The summation is composed of the two terms as ^2 q q ,^ — ^2 q i =k+q + Y^q'=ik+q • ^ is known that the dominant 
contribution arises when a particle interacts with itself and it belongs to the term q' = k + q. Then, we can separate 
the dominant contribution into two parts q = q' and q ^ q' again, but the term for q / q' vanishes using the 
argument of random phase approximation. Therefore, we have 



(0| 7)fe 7_ fc |0) ~ (0| 4 +9 aX- fe+9 'VlO> 

q' = k+q 

= (0| a q ,a q a q a q i |0) + (k— dependent minor terms). 

q=q' #o 

Within this approximation, using the n k in Eq. (fT5)) , we can write Eq. (|2"Tj) as 



(0| 7 fc7-fc|0) * £\ 



q^O 

E 

q^o 
N 



1 ( e„ + ng 



[k— dependent minor terms) 



V + 0(r, 2 ) 



(21) 



(22) 



where r) = ^na? /it <C 1. We put h = m = 1, and J2 q — (N/2n 2 n) J q 2 dq for convenience. The cutoff of the 
k-dependent minor terms is at most order of r] 2 or smaller. The effective range of the cutoff is judged the by an 
iteration method self-consistently. 

Substituting (0|7fc7_fc|0) in Eq. ([22|) into Eq. lfl9|) . we obtain a new higher order n k as 



n k = - 



1 ( e k + ng 



1 



8^ n 2 g 2 e k 2 
Et 



n + 0(n 2 ). 



(23) 



Now, let us back the new n k into Eq. ([22]) to obtain a new higher order (0|7fe7_fc|0), too. 



<0| 7 /c7- fc |0) = Y, 

qj&O 

= N 



ng 



E„ 



- 1 



7T 

3 



8\ n 2 g 2 e q ' 



V + 0(r, 2 ) 



V + 2 7T 



10 

y 



7r ) r) + 



16 



r] 3 + 0{n 3 ) 



(24) 



Therefore, we can trust up to rj 2 term in (0|7/ c 7_/ c |0) with the cutoff. 

Finally, substituting Eq. (|24[) into Eq. (|19[) , wc obtain the next higher order two coefficients C2 and C3. The particle 
depletion from the zero momentum condensate, (N — Nq)/N, is 



N E Uk 2N E 



fc#0 



£fc + 



- 1 



AT ^ 



22 2 
n z g z e k 



N Ei 



8 



7T--U + 2 7T-- --IT 0(jf) 



10 



-7? + 2vr 7T - - U 2 + 4?r 7T - - — - 7T h? + £>(?? ) 



10 



3^ 



i<"" 3 > 1/2+2 H) 



10 

Y 



7T J (na 3 ) 3 / 2 + G((na 3 ) 2 ). 



(25) 



Therefore, finally, we obtain c 2 = 2(tt - 8/3) and c 3 = (4/VtF) (tt - 8/3) (10/3 - tt). 
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V. SUMMARY 

We expanded the Hamiltonian of a homogeneous and interacting Bose system up to the yHVo terms, and applied 
the canonical transformation twice to find the average value of the higher order terms. An iteration method was 
applied to make up the cutoff. Two higher order coefficients c-i and C3 were obtained analytically. 
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